Draft version January 12, 2010 

Preprint typeset using WT^^ style emulatcapj v. 11/10/09 



HYBRID 7DORADUS - JSCUTI PULSATORS: 
NEW INSIGHTS INTO THE PHYSICS OF THE OSCILLATIONS FROM Kepler OBSERVATIONS 

A. Grigahcene\ V. Antoci^, L. Balona'\ G. Catanzaro^, J. Daszynska-Daszkiewicz'', J. A. Guzik'', G. Handler^, 

G. HOUDEK^ D. W. KURTZ'^, M. MARCONI^ M. J. p. F. G. M0NTEIR0\ A. MOYa'\ v. RIPEPI^ J.-C. SUAREZ^", 

K. Uytterhoeven", W. J. Borucki^^, T. M. Brown^'', J. Christensen-Dalsgaard", R. L. Gilliland^^, 
J. M. Jenkins^**, H. Kjeldsen", D. Koch^^, S. Bernabei^', P. Bradley^*^, M. Breger^, M. Di Criscienzo^^, 
M.-A. DuPRET^", R. A. Garcia", A. Garcia Hernandez^", J. Jackiewicz^\ A. Kaiser^, H. Lehmann^^, 
S. Marti'n-Ruiz^", p. Mathias^^, J. Molenda-Zakowicz'\ J. M. Nemec^^, J. Nuspl^^, M. Papar6^°, M. RoTH^^ 

R. SzABO^', M.D. SuRAN^^, R. Ventura'' 

Draft version January 12, 2010 

ABSTRACT 

Observations of the pulsations of stars can be used to infer their interior structure and test theoretical 
models. The main sequence 7Doradus (Dor) and JScuti (Set) stars with masses 1.2-2.5 Mq are 
particularly useful for these studies. The 7 Dor stars pulsate in high-order g modes with periods of 
order 1 day, driven by convective blocking at the base of their envelope convection zone. The 5 Set 
stars pulsate in low-order g and p modes with periods of order 2 hours, driven by the n mechanism 
operating in the He 11 ionization zone. Theory predicts an overlap region in the Hertzsprung-Russell 
diagram between instability regions, where 'hybrid' stars pulsating in both types of modes should 
exist. The two types of modes with properties governed by different portions of the stellar interior 
provide complementary model constraints. Among the known 7 Dor and 5 Set stars, only four have 
been confirmed as hybrids. Now, analysis of combined Quarter and Quarter 1 Kepler data for 
hundreds of variable stars shows that the frequency spectra are so rich that there are practically 
no pure 5 Set or 7 Dor pulsators, i.e. essentially all of the stars show frequencies in both the 5 Set 
and 7 Dor frequency range. A new observational classification scheme is proposed that takes into 
account the amplitude as well as the frequency, and is applied to categorize 234 stars as 5 Set, 7 Dor, 
(5Sct/7Dor or 7Dor/(5Sct hybrids. 

Subject headings: space vehicle - stars: oscillations - stars: variables: 5 Scuti - stars: variables: 
7 Doradus 
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1. INTRODUCTION 

Many stars oscillate in multiple simultaneous modes 
during some portion of their evolution. The frequency 
range, spacing, amplitude, or other properties can be 
used to infer the interior structure of these stars and test 
the input physics and predictions of theoretical models, 
a process called asteroseismology. The 7 Doradus (Dor) 
and S Scuti (Set) stars that pulsate in many simultaneous 
frequencies are particularly useful for asteroseismic stud- 
ies. These stars are core hydrogen-burning and some- 
what more massive than the Sun, at 1.2-2.5 Mq. How- 
ever, unlike the Sun, they have convective cores, shal- 
lower convective envelopes, and often rapid rotation. 

While heliose ismologists have be en searching for solar 
g modes (e.g., iGarcia et al.l 120071 ) . stellar astronomers 
have found them in several groups of pulsating stars in- 
cluding the 7 Dor variables. The 7 Dor g modes have 
periods of the order of one day, making them difficult to 
observe continuously from the ground. In addition, the 
g modes reach the stellar surface with small amplitudes 
which make their detection a challenge. 

From a pulsational point of view, there is a clear dis- 
tinction between S Set and 7 Dor stars. The former are 
short-period pulsating stars, with periods between 0.014 
and 0.333 day (d), while the latter are long-period pul- 
sators, with periods between about 0.3 and 3d. The 
distinction i s even clearer if one conside rs the pulsation 
constant Q (iHandler fc Shobbrookll200l : the <5Sct stars 
have Q < 0.055 d whereas for 7 Dor stars Q > 0.24 d. 
When the 7 Dor variables were first recogniz ed as a new 
group (|Balona et al. lll994l:lKave et al.lll999D . thefr posi- 
tion on the Hcrtzsprung-Russell diagram (HRD) already 
suggested some relationship with the 6 Set group. The 
7 Dor stars lie in a zone close to the cool border of the 
classical instability strip, partially overlapping with the 
(5 Set pulsators. This latter property led astronomers to 
investigate the possible existenc e of hybri d stars exhibit- 
ing b oth types of pulsations (jBreger fc Beichbuchncrl 
flOM ). 

These hybrid objects are of great interest because they 
offer additional constraints on stellar structure. The 
7 Dor stars pulsate in g modes, giving us the opportunity 
to probe the stellar core, while the S Set stars pulsate in 
p modes, that help to probe the stellar envelope. More- 
over, since 7 Dor stars pulsate in high-order g modes, 
the asymptotic approximation can be used, where we 
can benefit from an ana lytical solution of the equations 
(jSmevers fc Movall2007[ ). Hence, the existence of hybrid 
stars provides a unique opport unity to e xtend this ad- 
vantage to S Set stars ({Handle r fc Shobbro ok 2002) . 

Furthermore, in the same region of the HRD where 
the 6 Set and 7 Dor instability strips overlap, solar-like 
oscillations are also predicted to occur for ^Sct stars 
(|Houdek et al.|[T999HSa^adi et al.|[200l . So far, no de- 
tections of sola r-like oscillations in such hot stars have 
been reported (jAntoci et al.l 120091 ) , which is probably 
due to the detection limits imposed by earthbound mea- 
surements. Stochastic excitation as an excitation mech- 
anism in 7 Dor star s has also been searched for (e.g. see 
iPereira et al.l [20071 ), but so far not confirmed. Kepler 
data are expected to reach the precision needed to ob- 
serve stochastically excited pulsation. A positive or a 
significant null detection would set strong constraints on 
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Fig. 1. — Color-M agnitude diagram of 5 Set (plus signs , 
IRodn'guez et al.l [2000l) . 7 Dor (open cireles, IHenrv et all [2007]) 
and stars that show both typ es of pulsation (star symbols , 
IHenrv fc bWc! '20051: IRowe et a.i]\200B : lUvtterhoeven et al.l [2003 : 
[Han dler 2009, Rowe , priv. comm.). The Zero-Age Main Sequence 
(dashcd-dotted line, [Crawford "1973) as w ell as the ob served bor- 
ders of the S S et (solid lines, jR odrfgucz & Bregerll2001l ) and 7 Dor 
(dashed lines, IHandler & Shobbrookl 120021 ) instability strips, are 
indicated. Lines on the left indicate blue edges while those on the 
right give the red edges of the instability strips. 

t he parametrization o f the convective envelopes. 

IPoretti et all (|2009f ). using data from the French-led 
CoRoT mission, show that over most of the periodogram 
of HD 292790 (V = 9.5), the amplitude of the noise level 
is 0.01 mmag for an observational time span of 57 d. 
In Kepler observations of the combined Quarter and 
Quarter 1 data of a 9.5 mag star (time span of 44 d) the 
amplitude of the noise level is 0.001 mmag. The noise 
level in Kepler observations is therefore approximately 
10 times smaller than in CoRoT observations. 

The large number of stars to be studied by Kepler 
will also allow for control of long-term (low-frequency) 
instrumental changes; understanding these is important 
for confident identification of the low- frequency g modes. 
Most importantly, the 3.5-y time span of the mission 
will give data with unprecedented frequency resolution. 
This is known to be necessary from ground-based studies 
that show closely spaced frequencies in data sets span- 
ning years. 

2. CURRENT OBSERVATIONAL STATUS 

Among the hundreds of known 5 Set stars and sev- 
eral dozen confirmed 7 Dor variables, only 4 have 
been previously observationally identified as hybrids 
and in terpreted afterward using theoretical models 
(|Bouab id et al. 2003). 

One of the major contributi ons came from th e Cana- 
dian Sp ace Mission MOST (jMatthewsl l2007f ). while 
CoRoT (jAuvergne et al.|[2009f ) is expected to add new 
confirmed cases in the near future, especially since 
the preliminary results indicate the possible detection 
of hybrid 7 Dor an d 5 Set pulsators in the exofield 
(jMathias et al.ll2009[ ). 

The known hybrid pulsators are depicted, together 
with single-type pulsators , in Fig. [J The first cas e 
recorded was HD 209295 (|Handler fcShobbrookl [200l) . 
but it turned out afterward that this star is a close 
binary and the long period modes could be tidally 
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Fig. 2. — Hybrid theoretical instability domain for £ = — 3 g 
and p modes. Models have: Z = 0.02, omlt = 2.0, aov = 0.2. 
The vertical line corresponds to the effective temperature of the 
Kepler target: KIC 11445913 = 7200 ± 200 K). 



excited ([Handler et al.l |2002[ ). Following this work, 
iHenrv fc Fekell ()2005[) ~discovered both 7 Dor a rid (5 Set- 
type pulsations in the single Am star HD 8801 |Handler| 

mm ). 

MOST fou nd two a d dition al hybrid pulsators: 
HD 114839 (iKing et al.l [20061 ) and BD+18 4914 
owe "eralll20060 . Both of these are also Am stars. In 
repa ration for the CoRoT mission. lUvtterhoeven et al.l 
2008( 1 identified HD 49434 as a candidate hybrid pul- 
sator. 

It is very interesting to note that three of the four 
known hybrids show Am type peculiarities, even though 
Am stars are less likely to show pul sational varia bility 
than A/F-type main-sequence stars (jKurtd [l989[ ). In- 
creasing the number of hybrids is expected to shed light 
on the relationship between hybrid pulsation and Am 
chemical peculiarity. 

3. THEORETICAL EXPECTATIONS AND MODELING 

It is well known that the K-mechanism operating in 
the Hen ionization zone is responsible for the pulsa- 
tional driving of (5 Set stars (| C he valiei] 119711 ). With this 
mechanism, the theoretical blue edge of the S Set instabil- 
ity strip is predicted correctly with linear non-adiabatic 
models for radial as well as for non-radial modes. But 
determining the theoretical red edge is more difficult, be- 
cause it requires a time-dependent treatment of the inter- 
action between convection and pulsation. In this stage 
of their evolution, cool S Set stars have relatively deep 
surface convection zones which typically extend beyond 
the location of the second stage of helium ionization, 
the principal driving region of 6 Set pulsations. Conse- 
quently, mode stability is crucially affected by convection 
dynamics for stars located near the red edge of the clas- 
sical instability strip in the HRD. The location of the 
red edge can be described theoretically only by means 
of a time-dependen t treatment of the turbulent fluxes 
(e.g.,|P cuprce"1977|: IBaker fc Gough1[T979t IGonzillMl 
IStelhn gwcrf 1984). Moreover, the observed location of 
the instability strip provides independent information for 
calibrating convection models. 

Several authors (e.g.. iBono et al.lll999l: iHoudekl [2OOOI: 
iDupret et all 120051 : IXiong fc DeiiJ I2007D have success- 
fully modeled the red edge of the classical instability strip 



but reported that different physical mechanisms are re- 
sponsible for stabilizing the pulsations. It is therefore of 
great importance that our calculations for mode stabil- 
ity be reassessed and Kepler is in the position to pro- 
vide the necessary high-quality data for calibrating and 
then testing the models for convection in classical pul- 
sators. Moreover, analysis of the follow-up ground-based 
multi-color photometry can yield valuable constraints on 
eflnciency of the convective transpo rt in these pulsators 
(jPaszvhska-Daszkic wicz et al.l 120031 ). 

[Guzi k et a l. (20 Q0|). us i ng fr ozen-convection models, 
and iDupret et ali (120051 l2006f). u sing time-dependent 
convection (|Grigahcene et al.ll2005D . have shown that the 
position of the convective envelope base is the key to driv- 
ing 7 Dor g modes. The convective blocking mechanism, 
operating when the convective timescale becomes com- 
parable to the g mode period, can explain the location 
of the blue edge of the 7 Dor instability strip and its sen- 
sitivity to the mixing length parameter umlt- On the 
other hand, the red edge of the 7 Dor instability strip is 
caused by radiative damping in the g mode cavity, which 
dominates over the excitation near the convective enve- 
lope base. As an example, we show in Fig. [2] the hybrid 
theoretical instability domain for I = 0—3 g andp modes. 
For the calculated grid (^=0.02, aMLT=2.0), the effec- 
tive temperature range of hybrids is 6900 K < Tcff < 
7500 K. The vertical line in Fig. [2] represents a model of 
a hybrid pulsator for which we have unstable high order 
g modes, corresponding to 7 Dor range as well as unsta- 
ble low order p and g modes, corresponding to typical 
5 Set pulsations separated by a gap of stable modes. 

We also note that the stability of modes in the gap 
decreases as degree ^ increases. This is due to the fact 
that the Lamb frequency Sf oc i{t -1-1), so the size of 
the propagation cavity, of g modes increases with i for a 
given frequency. 

Therefore, the direct comparison of the observed fre- 
quency spectrum with models is a test of the excitation 
models and can improve our knowledge about important 
physical characteristics of the stars, including the age, 
chemical composition, etc. 

Another more general problem in stellar pulsation 
modeling is connected with the effects of rotation on 
eigenfrequencies and eigenfunctions. The commonly 
adopted approach is perturbation theory which assumes 
that oscillation frequencies are much larger than the an- 
gular rotation rate VL and that the star is not signifi- 
cantly deformed by the centrifugal force (O <C ilcrit = 
y/GM]W) (jReese et al.ll2008f ). This assumption is vi- 
olated for many 5 Set and 7 Dor stars, due to their rel- 
atively rapid rotation. This implies that rotation will 
significantly affect the pulsation frequency spectrum and 
surface amplitude distribution. This in turn will affect 
how limb darkening and light cancellation affects the pul- 
sation mode visibility and radial velocity variations. 

Hybrids can be extremely useful when investigating the 
internal rotation profile of stars. In[Suarez et al. (200^) 
the authors examine the effect of rotation on the oscilla- 
tion spectrum for intermediate-mass rotating stars. Sig- 
nificant effects for both p and mixed modes (very low 
radial order g and p modes) arc predicted. Analysis of 
hybrid stars will allow us to constrain the modeling of 
their rotation profiles. Better constraints on the internal 
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rotation profile will help us better understand angular 
momentum transport processes in the stellar interior. 

Neither 5 Set stars nor hybrid stars are expected to 
pulsate in asymptotic p mode o s cillati ons. However, 
recently IGarcia-Hernandez et al.l (|2009D found regular 
spacings in the oscillation spectrum of the CoRoT S Set 
star HD174936. Analysis of the 422 detected osciUa- 
tion frequencies revealed a spacing periodicity of around 
52 /zHz. Although the modes considered were not in 
the asymptotic regime, a comparison with stellar models 
confirmed that this signature may stem from a quasi- 
periodic pattern similar to the so-called large separa- 
tion in solar-like stars. Other frequency spacings have 
also been obser ved in other stars using ground based ob - 
servations (e.g. [Handler et a"I1[2000l: iBreger et al.l [20091 ). 
From numerical results, this approximate asymptotic be- 
havi our has been found to hold fo r p modes of low order 
(e.g. lChristensen-Dalsgaardll20do| ) . The presence of such 
spacings in hybrid stars may allow us to analyze simul- 
taneously p and g modes with some characteristics of 
the asymptotic regime and even solar-like oscillations, if 
detected. 

Interestingly, asteroseismology of hybrid stars can also 
benefit from techniques used for analyzing and interpret- 
ing the oscillation spectra of both type of pulsations. For 
instance, if g modes in the asymptotic regime are de- 
tected, techniques based on the asymptotic properties of 
g modes can be applied, namely t he Frequency Rati o 
Method (FRM - IMova et all [20051 : IS^ez et al.l[2005h . 
The FRM is particularly adapted for obtaining astero- 
seismic information on 7 Dor pulsating stars showing at 
least three pulsation frequencies. The method provides 
an identification of the radial order n and degree £ of 
observed frequencies and an estimate of the integral of 
the buoyancy frequency (Brunt- Vaisala) weighted over 
the st ellar radius along the radiative zone. iMiglio et all 
()2008D also considered high-order g modes to probe the 
properties of convective cores in Slowly Pulsating B and 
7 Dor stars. 

4. KEPLER'S FIRST OBSERVATIONS 

The main result from the inspection of this first Ke- 
pler data is that there are practically no pure 5 Set or 
7 Dor pulsators. This finding raises a new problem: why 
is it that there is a clear distinction between the two 
types of pulsator in ground-based observations but not 
in Kepler observations? Part of the answer could be that 
in ground-based photometric observations the modes of 
highest visibility are modes of low degree whereas at the 
photometric precision of Kepler many more modes of 
high degree are now visible. 

Now we must consider whether we can construct a 
classification scheme that allows us to retain the useful 
5Sct/7Dor categories. It is important to notice that re- 
lying only on the mode frequencies, Kepler^s data show 
that all stars with variability in these frequency ranges 
are hybrids. The use of amplitude in addition to fre- 
quency can help by potentially filtering out the high- 
degree modes that generally would be expected to have 
a lower photometric amplitude. Considering both am- 
plitude and frequency, we propose a new observational 
classification scheme for S Set and 7 Dor stars: 

• S Set: most of the frequencies are > 5 d^^, and the 
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Fig. 3. — Theoretical HRD showing classification of Kepler tar- 
get stars. Filled circles are (5 Set; open circles S Sct/jDoi; crosses 
7 Dor /(5 Set and plus signs 7 Dor. The solid lines show the Zero- 
Age Main Sequence and the radial fundamental red and blue edges 
IR, IB) and the 4- th overtone radial red and blue edges {4R, 4B) 
Dupret et al.|[2005l) . The dashed lines are the red and blue edges 
of th e 7 Dor instability strip {£ = I and mixing length auLT = 
2.0) HDupret et al.ll20n5l) . 



TABLE 1 

Number of stars and percentage in each glass. For each 
class the mean values of the effective temperature is 

GIVEN. 



Class 


Number 


Percent 


< log Teff > 


<5Sct 


67 


27 


3.885 ± 0.003 


(5Sct/7 Dor 


32 


14 


3.883 ± 0.006 


7Dor/<5 Set 


19 


9 


3.868 ± 0.006 


7 Dor 


116 


50 


3.853 ± 0.005 



lower frequencies are of relatively low amplitude; 

• i5Sct/7Dor hybrid: most of the frequencies are > 
5 d~^, but there are some lower frequencies which 
are of comparable amplitude; 

• 7 Dor: most of the frequencies are < 5 d~^, and the 
higher frequencies are of relatively low amplitude. 

• 7 Dor/ J Set hybrid: most of the frequencies are < 
5 d~i^, but there are some higher frequencies which 
are of comparable amplitude; 

We apply this scheme to a sample of 554 stars selected 
from the Kepler working group targets, where we have 
excluded stars that are constant or eclipsing variables. 
Among them, 234 show 7 Dor or 5 Set frequencies and 
have an effective temperature near the instability strips 
for these types of pulsation. Table 1 summarizes the 
categorization and mean properties for these stars. The 
same stars are shown in the theoretical HRD in Fig. m 
using the Toff and radii derived from the Kepler Input 
Catalogue (KIC). Also shown a re calculated red an d blue 
edges of the instability strips (jDupret et al.l 120051) . We 
see that there is some separation between categories, but 
it is not as clean as one might have hoped. We also see 
that the theoretical red and blue edges are not a good 
match for the observations. 



Hybrid 7 Dor ~ d Set pulsators from Kepler 



5 



Frequency [|.iHz] 




9 12 15 

Frequency [d"'] 



Frequency [nHz] 

120 150 180 210 240 270 




12 

Frequency [6''] 



Fig. 4. — Frequency spectra of the long-cadence data obtained 
by Kepler for the stars KIC 9775454 and KIC 11445913 that dis- 
play hybrid pulsations in the 7 Dor region (below 5d~^) and in the 
(5 Set region (above 12 d~^). The details of the spectra at higher 
frequencies require short-cadence data to be confirmed and ana- 
lyzed. 

The amplitude spectrum of two hybrid stars 
(KIC 9775454 and KIC 11445913) are shown in Fig. H 
The vertical line in Fig. [2] indicates the position in effec- 
tive temperature of the star KIC 11445913, confirming 
it as a potential hybrid candidate. Fig. 2] also shows the 
very low level of noise measured in the gap between the 
frequency ranges for the p modes and g modes. This 
gap, corresponding to a lack of unstable modes, is ex- 
pected in the theory for hybrid pulsators. However, in 
some stars observed by CoRoT, several frequencies can 
be found in the "gap" . The same behavior is also found 
for some of the Kepler targets. An explanation for these 
frequencies still needs to be found from a more detailed 
study of the observations and by modeling the stars be- 
ing observed. In some cases, rotational splitting seems 
to be the cause of modes in the freq uency gap between 
p and g modes (jBouabid et al.|[2009t) . Preliminary anal- 
ysis of Kepler data for candidate hybrid stars gives very 
encouraging results for the number of hybrid stars with 
very high quality frequency spectra. The large number of 
pulsational frequencies makes it possible to study global 
characteristics of the frequency spectrum such as asymp- 
totic properties, and also individual frequency fittings. 

5. CONCLUSIONS 

Our first look at the Kepler data provides encouraging 
results on the number of hybrid stars thanks to very high 



quality frequency spectra coupled with unprecedented 
photometric precision. The Kepler data are essential 
to overcome the aliasing and photometric stability prob- 
lems affecting the 7 Dor frequency range associated with 
ground-based photometric campaigns. 

The additional frequencies found in hybrid stars will 
help to overcome difficulties with mode identification and 
the detection of many 7 Dor g modes will enable analy- 
sis by asymptotic techniques. The Kepler data on 7 Dor 
and S Set stars promise to help us to answer a number of 
outstanding questions: why have Kepler (and CoRoT) 
detected a significant number of hybrids, whereas theory 
predicts the existence of hybrids in only a small over- 
lapping region of the instability strips? What is the 
connection between Am chemical peculiarity and hybrid 
behavior? Why are modes observed between the S Set 
and 7 Dor frequency ranges, where theory predicts a fre- 
quency gap? Will we be able to confirm stochastic ex- 
citation or other pulsation driving mechanisms? Can 
the high-photometric precision of Kepler detect higher- 
degree modes {i>2) that could fill in the frequency gap? 
Will these stars help us to distinguish between explana- 
tions for the red edge of the S Set instability strip? Will 
we be able to discern rotational frequency splittings and 
derive internal differential rotation rates for these stars? 

Perhaps just as important are the many Kepler targets 
that lie within the 7 Dor or i5 Set instability strips that 
show no periodic variations at micromagnitude photo- 
metric precision levels obtainable by Kepler. Monitoring 
over several years with Kepler will provide insight on 
amplitude limiting and mode selection mechanisms for 
these stars that is not possible with ground-based obser- 
vations. 

The analysis presented here of the first and second 
Quarter data obtained by Kepler confirms its capabil- 
ity to provide very high quality data that will allow us 
to study hybrid pulsators of the S Set /7 Dor type. More- 
over, the very rich frequency spectra of these stars leads 
us to propose a new classification scheme to overcome 
the ambiguity of the older, pure frequency-based scheme. 
The close frequency spacing of relatively short period 
modes in 6 Set stars requires short-cadence data. We ex- 
pect to obtain short-cadence data for selected targets of 
mixed character in the near future. 
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